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Abstract:  -  For  the  semilinear  elliptic  equation  Au  +  p(x)u~'*  =  0,  x  6  Rn, 
n  >  3,  7  >  0,  we  show  via  the  barrier  method  the  existence  of  a  positive 
entire  solution  behaving  like  |a'|2_n  near  co. 

1      Introduction 

We  study  the  singular  semilinear  elliptic  equation 

(1)  Autp(x')w_1  =  0 

in  Rn.  This  type  of  equation  arises  in  the  boundary  layer  theory  of  viscous 
fluids  [3,4].  From  the  results  of  Fulks  and  Maybee  [8],  Crandall,  Rabinowith, 
and  Tartar  [5],  Gomes  [9],  and  recently  Lazer  and  McKenna  [14],  it  follows 
that  (1)  has  a  unique  classical  solution  within  a  bounded  domain  Q,  where 
p(x)  is  a  sufficiently  regular  function  which  is  positive  on  Q.  Kusano  and 
Swanson  [12]  gave  the  existence  proof  on  exterior  domains.  As  for  the  ex- 
istence of  entire  solutions,  not  much  is  known.  Edelson  [7],  Kusano  and 
Swanson  [13]  have  been  able  to  show  the  existence  of  entire  solutions  of  (1) 
with  7  E  (0, 1),  and  p(x)  sufficiently  regular.  In  this  paper  we  show  via  the 
upper  and  lower  solution  method,  which  is  also  referred  to  as  the  barrier 
method,  that  (1)  has  a  bounded  positive  entire  solution  vanishing  at  oo  in 
Rn  for  n  >  3  and  all  7  >  0. 

The  author  learned  after  this  paper  was  finished  that  a  similar  result  was 
given  earlier  by  R.  Dalmasso  [6],  but  by  a  different  approach. 
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2     Preliminaries 

We  first  state  the  theorem  by  Kusano  and  Swanson  [13]  for  the  case  0  <  7  < 
1. 

Lemma  1.  Equation  (1)  has  an  entire  bounded  positive  solution  u(x)  in  Rn 

for  n  >  3  ,and  \x\n~2u(x)  is  bounded  and  bounded  away  from  zero  near  00  if 

p(x)  satisfies  the  following  conditions: 

(HI)  p(x)  £  C?0C(Rn),  n  >  3,  p(x)  >  0,  x  £  JT  \  {0}, 

(H2)  3C  >  0,  such  that  C<f>(\x\)  <  p(x)  <  <j>{\x\),  <f(x)  =  max\x\=tp(x),  0  < 

t  <  00, 

(HS)  /r  tn-1+^n-V<f>{t)dt  <  00. 

The  term  "entire"  has  often  been  used  for  solutions  of  equation  (1)  in 
Rn.  To  avoid  confusion  with  the  traditional  definition  for  entire  functions, 
we  use  the  term  "C2+Q-entire".  A  C2+Q-entire  solution  of  (1)  is  defined  to 
be  a  function  u(x)  G  C?£.a(Rn)  that  satisfies  (1)  pointwise  in  Rn. 

The  method  that  we  shall  be  using  heavily  in  our  proof  is  the  so-called 
barrier  method,  or  upper-lower  solution  method. 

We  consider  the  elliptic  boundary  value  problem 

/«\  (  Lu  +  f(x,u)  =  0  in  D 

^  '  1    Bu  =  adujdu  +  bu  =  g    on  3D 

where  D  is  a  smoothly  bounded  domain  in  RN  and  v  =  (^i,---,^n)  is  a 
smoothly  varying  outward  normal  vector  field  on  dD  which  is  of  class  C2+a, 
while  a  and  b  are  positive  constants.  We  also  assume  that  /  £  Ca  and  that 
g  has  an  extension  g  to  the  interior  of  D  such  that  g  £  C2+Q . 

An  upper  solution  to  the  above  problem  is  a  function  <f>  satisfying 

f  L<f>  +  f{x,<f>)  <  0      in  D 
\  B4>>g  on  dD. 

A  lower  solution  to  the  above  problem  is  a  function  ip  satisfying 

f  Lip  +  f(x,il>)  >  0      in  D 
\  B^<g  on  dD. 

We  assume  that  dD,f,g,  and  the  coefficients  of  L  are  smooth  in  what 
follows. 


Lemma  2.  (Theorem  2.3.1  of  [16])  Let  (f>  be  an  upper  solution  and  t/>  a 
lower  solution  with  i\>  <  (f)  on  D .  Then  there  exists  a  solution  u  to  the  above 
boundary  value  problem  with  tp  <  u  <  <f>. 

We  consider  the  following  example: 

j  u"  +  Xu  -  u3  =  0    x  £  (0,tt) 

1  U  —  0  X  =  0,  7T 

By  the  above  theorem,  if  A  >  1,  then  the  problem  has  at  least  three  solutions. 

Actually,  u  =  esin  x  with  e  small  is  a  lower  solution,  and  u  =  Rxll2  with 
R  large  is  an  upper  solution.  Therefore  there  exists  a  solution  u  such  that 
u  <  u  <  u  in  (0,7r).  Clearly  —  u  and  0  are  also  solutions  to  this  problem. 

The  following  lemma  on  the  barrier  method  for  D  —  Rn  is  due  to  Ni 
[15]  in  1982.  A  special  case  was  proved  earlier  by  Ako  and  Kusano  [1]  in 
1964.  The  proof  is  standard.  Using  the  well  known  result  on  the  upper-lower 
solution  approach  in  bounded  regions  (see  Sattinger  [16]),  we  first  solve  the 
equation 

Lu  +  F(x,u)  =  0 

on  Br.  Then  by  letting,  R  — *  oo,  we  obtain  a  solution  on  Rn  by  a  diagonal 
process. 

Lemma  3.  Let  ux  >  u2  in  Rn  be  such  that 

(  {  Iwi  +  f(x,ui)  <  0 

1   '  1  Lu2  +  f(x,u2)>0 

where  f  is  locally  Holder  continuous  in  (x,u)  and  locally  Lipschitz  in  u,  and 
L  is  an  elliptic  operator  of  second  order.  Then  there  exists  a  solution  u  of 
Lu  +  /(#, u)  —  0  with  Ui  >  u  >  u2. 

3     Main  Result 

Theorem  1.  Under  the  same  conditions  as  given  in  Lemma  1,  the  equation 
(1)  has  a  C2+a-entire  positive  solution  in  RN ,N  >  3,  vanishing  at  oo  at  the 
rate  of  at  least  Ixl9^-2^  with  some  q  £  (0, 1)  for  any  7  >  0. 


The  difficulty  in  constructing  the  proof  is  to  find  an  appropriate  upper 
solution  for  equation  (1).  In  order  to  use  the  barrier  method  we  first  study 
the  nonsingular  equation 

Au  +  p(x)[6  +  u]~7  =  0. 

For  each  fixed  7  there  corresponds  a  solution  u7(rr).    Letting  7  — *  00,  we 
show  that  the  limiting  function  is  the  desired  solution. 
Proof:    By  Lemma  1,  for  7  =  71  G  (0,1),  equation  (1)  has  a  C2+a-entire 
positive  solution  ux(x)  in  i?n,  n  >  3,  vanishing  at  00  at  the  rate  r2_n.   We 
claim  that  u  =  cu\  is  an  upper  solution  of  the  equation  (1)  for  7  >  1,  where 

q  <-r— —  <  1, 
1  +7 

^1+71-9(1+7)     _j_ 

c>( )1+"S    M  =  maxx€Rn\u(x)\. 

Q 

In  fact: 

AS  +  ^1 

=  cq{q  -  l)uq-2\X7u\2  -  cqug~l p(x)u~11  +  pix^u'^ 
<  —  cquq~1p(x)u~11  -f  p{x)c~1u~iq 

_      PJX)    ^    _  ^^Q  N 

C7W79^  Ul+7l -9(1+7)  ^ 

<  p(x)       _     M^-^+^)  q 

—  ciui  q  ^1+71-9(1+7)' 

<pW(1_1)  =  0. 

Let  6  be  a  fixed  positive  number.  We  then  observe  that  u  is  an  upper 
solution  of  the  equation 

(4)  Au{x)  +  p{x)[u{x)  +  S}-1  =  0,  x  G  iT. 

u  =  0  is  a  lower  solution  of  (4).  Since  u  =  cu\  >  0,  u  >  u  in  Rn.  By  Lemma 
2,  (4)  has  a  solution  u  such  that  u  <  u  <  u. 


For  6  <  6,  u  is  a  lower  solution  of  (4)  with  6  =  6.  Lemma  2  then  implies 
that  (4)  has  a  solution  u  for  6  =  6  such  that  u  <  u  <  u. 

Let  {6n}f  be  a  sequence  of  strictly  decreasing  positive  numbers,  and  let 
un(x)  be  a  smooth  positive  solution  of  (4)  when  6  =  6n.  From  the  construc- 
tion of  our  lower  solutions,  it  is  clear  that  un(x)  >  un_1(x)  for  all  n.  So 
limn^oo  un(x)  =  u(x)  exists  for  all  x  G  Rn  and 

(5)  u  <  u  <  u 

for  x  E  Rn. 

We  can  now  assert  that  u  G  C2+a(Rn)  and  that 

(6)  Au  +  p(x)u"r  =  0 

for  x  E  Rn.  This  follows  from  more  or  less  standard  arguments. 

Let  x0  E  Rn  and  r  >  0.  We  consider  the  ball  of  radius  r  centered  at  x0, 
B(x0,r)  in  i?n.  Let  $  be  a  C00  function  which  is  equal  to  1  on  B(x0,r/2) 
and  equal  to  0  off  B(x0,r).  We  have 

A(*un)  =  2W-Vun+pn 

for  77  >  1,  where  pn  is  a  term  whose  L°°  norm  is  bounded  independently  of 
n.  Therefore  for  n  >  1  we  have 

where  6nj-,j  =  l,---,n  and  qn  are  terms  bounded  independently  of  n  for 
n  >  1.  Integrating  the  above  equation,  we  have  that  there  exist  constants 
C\  >  0  and  c-i  >  0  independent  of  n  such  that 

/  |Vun|2^<Cl(/  \Vun\2dx)^2  +  c2. 

JB{xo,r)  JB(x0,t) 

From  this,  it  follows  that  the  L2(B(x0:  r))-norm  of  |V$itn|  is  bounded  inde- 
pendently of  n.  Hence,  the  L2(B(x0,r/2))-norm  of  |Vun|  is  bounded  inde- 
pendently  of  n.  Let  tyy  be  a  C°°  function  which  is  equal  to  1  on  B(x0,r/4) 
and  equal  to  0  iff  B(x0,r/2).  We  have  for  n  >  1, 

A(*1un)  =  2V*1.V«„+plnj 


where  pln  is  a  term  whose  L°°(B(x0,r/2))-  norm  is  bounded  independently 
of  n.  From  standard  elliptic  theory,  the  W2,2(i?(:E0,r/2))-norm  of  \&iun  is 
also  bounded  independently  of  n  and  hence,  the  W2'2(B(x0,r/4))-noTm  of 
un  is  bounded  independently  of  n.  Since  the  Wl'2(B(x0,r /4))-norm  of  the 
components  of  Vun  are  bounded  independently  of  n,  it  follows  from  the 
Sobolev  embedding  theorem  that  if  q  =  2n/(n  —  2)  >  2  for  n  >  2  and  in 
addition  if  q  >  2  is  arbitrary  for  n  <  2,  then  the  L9(5(x0,r/4))-norm  of  |wn| 
is  bounded  independently  of  n.  Let  \&2  be  a  C°°  function  which  is  equal  to  1 
on  B(x0,r/8)  and  equal  to  0  iff  B(x0,r/4).  We  have  for  n  >  1, 

A(tf2un)  =  2Vtf2-Vtin+p2n, 

where  p2n  is  a  term  whose  L°°(B(x0,  r/4))-  norm  is  bounded  independently  of 
n.  Since  the  right  hand  side  of  the  above  equation  is  bounded  in  Lq(B(x0,r/4)) 
independently  of  n,  the  W2'9(B(x0,  r/4))-norm  of  $2un  is  also  bounded  inde- 
pendently of  n.  Hence,  the  W2'q(B(x0,r /8))-norm  of  un  is  bounded  indepen- 
dently of  n.  Continuing  this  line  of  reasoning,  after  a  finite  number  of  steps, 
we  find  a  number  ra  >  0  and  q\  >  n/(l  —  a)  such  that  the  H/2'9l(5(x0,r1))- 
norm  of  un  is  bounded  independently  of  n.  Hence,  there  is  a  subsequence 
of  {wn}i°'  which  we  may  assume  is  the  sequence  itself,  which  converges  in 
C1+a(B(a;0,r1)).  If  6  is  a  C°°  function  which  is  equal  to  1  on  (JB(rr0,r1/2)) 
and  0  off  B(x0:r1),  then 

A(9un)  =  2V0  •  S7un  +  pn,  where  pn  =  6Aun  +  unA6. 


The  right-hand  side  of  the  above  equation  converges  in  CQ(B(x0,  ri)).  Hence 
by  Schauder  theory,  {9un}^°  converges  in  C2+a(B(x0,ri))  and  thus  {un}i° 
converges  in  C2+a(B(x0,ri/2)).  Since  x0  was  arbitrary,  this  shows  that  u  E 
C2+a(Rn).  Clearly  (6)  holds.  ' 

4     Some  remarks 

Remark  1: 

For  n  =  1,  the  properties  of  positive  solutions  of  equation  (1)  have  been 
studied  by  Taliaferro  [17],  and  Gatica  [10].    For  n  =  2,  no  entire  positive 
solution  of  equation  (1)  exists  regardless  of  its  asymptotic  behavior  at  oo 
(see  [13]). 
Remark  2: 


It  is  observed  by  Callegari,  Friedman  and  Nachman  [2],  [3,4]  that  if  the 
partial  differential  equations  describing  the  boundary  layer  behind  a  rarefac- 
tion or  shock  wave  (with  viscosity  proportional  to  the  temperature)  traveling 
down,  and  perpendicular  to,  a  flat  plate  are  written  in  terms  of  a  stream  func- 
tion and  a  similarity  variable  the  following  Blasius-type  equation  emerges 
[18]. 

/'"M  +  /WfM  =  o, 

where 

/(0)  =  0,    f'(Q)  =  K,    /'(oc)  =  l. 

Here,  0  <  K  <  1,  for  rarefaction  waves  and,  1  <  K  <  6,  for  shock  waves. 
(K  =  0  corresponds  to  the  classical  Blasius  problem.)  Adopting  the  Crocco 
variables 

results  in  the  system 

gg"  +  x  =  0, 

g'(K)  =  0,   g(l)  =  0, 
which  falls  into  the  class  of  equation  discussed  in  this  paper. 
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